Complete hierarchies of efficient approximations to problems in entanglement theory 
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We investigate several problems in entanglement theory from the perspective of convex optimization. This 
list of problems comprises (A) the decision whether a state is multi-party entangled, (B) the minimization of 
expectation values of entanglement witnesses with respect to pure product states, (C) the closely related evalua- 
tion of the geometric measure of entanglement to quantify pure multi-party entanglement, (D) the test whether 
states are multi-party entangled on the basis of witnesses based on second moments and on the basis of linear 
entropic criteria, and (E) the evaluation of instances of maximal output purities of quantum channels. We show 
that these problems can be formulated as certain optimization problems: as polynomially constrained problems 
employing polynomials of degree three or less. We then apply very recently established known methods from 
the theory of semi-definite relaxations to the formulated optimization problems. By this construction we arrive 
at a hierarchy of efficiently solvable approximations to the solution, approximating the exact solution as closely 
as desired, in a way that is asymptotically complete. For example, this results in a hierarchy of novel, efficiently 
decidable sufficient criteria for multi-particle entanglement, such that every entangled state will necessarily be 
detected in some step of the hierarchy. Finally, we present numerical examples to demonstrate the practical 
accessibility of this approach. 
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I. INTRODUCTION 

One of the reasons for the superior performance of quan- 
tum devices for computation and communication compared 
to their classical counterparts is simply due to the fact that 
in quantum mechanics, one has a very large space at hand 
to work with: the dimension of the state space of a number of 
quantum bits is expontially larger than the corresponding con- 
figuration space of classical bits. This renders the simulation 
of a quantum computer on a classical device a difficult task. 
But it is not only the sheer size of state space that makes the 
assessment of quantum states a difficult problem. In fact, even 
to decide whether quantum states have certain properties that 
are of central interest in quantum information science often 
amounts to solving computationally hard problems on a clas- 
sical computer. Most prominently, to decide whether a known 
state p of a finite-dimensional bi-partite system is separable 
or entangled, i.e., whether or not it can be written as a convex 
combination of product states 

n 

is already an NP-hard problem in the system size [jj]. A state 
is separable if there is a preparation of the state that involves 
only local quantum operations and shared classical random- 
ness. Such states are correlated, but classically correlated, as 
the source for the correlations can be thought of as resulting 
entirely from the shared source of randomness 0]. Due to 
the central status of the concept of entanglement in quantum 
information, a very significant amount of research has been 
dedicated to the problem of finding good criteria for separa- 



bility that are suitable for specific contexts 01 ■ 

To state whether a state is separable or not is equivalent to 
stating whether a state is in the convex hull of product states. 
Also, the evaluation of many measures of entanglement es- 
sentially require the solution of a convex problem. So in re- 
cent years, it has increasingly been realized that a good deal of 
new insight in several problems in quantum information and in 
particular in entanglement theory could in fact come from the 
field of research that is primarily concerned with questions of 
this type fllHIMIcll/TMM[loll : this is the theory of convex op- 
timization. Many problems are already of the required form, 
and powerful tools such as the concept of Lagrange duality 
readily deliver bounds for the problems at hand. Examples 
include the evaluation of measures of entanglement that rea- 
sonably quantify the degree of entanglement of a given state, 
such as the distillable entanglement or the asymptotic relative 
entropy of entanglement 0y|]. Also, it has been realized that 
while the complete solution of the question of separability is 
NP-hard, one can nevertheless find hierarchies of sufficient 
criteria for entanglement in the bi-partite setting. In each step, 
by solving an efficiently solvable convex optimization prob- 
lem, one finds an answer to the problem in the form (i) one 
can assert that the state is entangled, or (ii) one cannot assert 
it, and has to go one (computationally more expensive) step 
further (||. The problem of testing for multi -partite entangle- 
ment has been related to robust semi-definite programming 
and a hierarchy of relaxations in Ref. |9[. 

This paper is concerned with a link of the theory of en- 
tanglement to the theory of convex optimization in a similar 
spirit. The central observation of this paper is very simple yet 
potentially very useful: many problems related to entangle- 
ment can be cast into the form of optimization problems with 
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polynomial constraints of degree three. This includes the (A) 
question whether a state is entangled or not, notably not only 
in the bi-partite, but also for the several separability classes of 
the multi-partite setting. Then, (B) the construction of non- 
decomposable witnesses involves a problem of this kind, as 
well as the (C) evaluation of the geometric measure of en- 
tanglement to quantify multi-partite entanglement. (D) Also, 
when considering entanglement witnesses based on second 
moments rather than on first moments one has to solve a prob- 
lem of this form. We will also discuss criteria based on linear 
entropies (i.e., p-norms forp = 2). (E) Finally, we will briefly 
mention the evaluation of maximal output purities of quantum 
channels with respect to p-norms for p = 2. This structure is 
due to the fact that in all these instances, one essentially min- 
imizes over product state vectors of a multi-partite quantum 
system. 

This polynomial part of the optimization problems is still 
non-convex and computationally expensive to solve. Yet, ap- 
plying r esults from relaxation theory of non-convex problems 
notably the method of Lasserre ll3ll . we 
find hierarchies of solutions to our original problems, and each 
step is a better approximation than the previous one. Each step 
itself amounts to solving an efficiently implementable semi- 
definite program fl^l . Moreover, the hierarchy is asymptot- 
ically complete, in the sense that the exact solution is asym- 
potically attained. The increase of the size of the vector of 
objective variables of these semi-definite problems grows no- 
tably polynomially in the label of the hierarchy. 

We will first clearly state how one can introduce auxiliary 
variables to cast the considered problems from entanglement 
theory into the desired form. Then, we will investigate the 
hierarchies of relaxations in detail, and study numerical ex- 
amples. Finally, we will summize what has been achieved. 



II. PROBLEMS IN ENTANGLEMENT THEORY AS 
OPTIMIZATION PROBLEMS 

The problems that we will encounter are of the following 
type or similar. At the core are typically minimizations over 
product vectors, originating from the very definition of the 
concept of entanglement. Given a W = W\ we seek the 
minimum of 

tr[|^i><Vi|®...® \^n){^n\W], (2) 

where the minimum is taken with respect to product state 
vectors of a composite quantum systems with parts labeled 
1 , . . . , N, with Hilbert space Ji = Ti\ ® .. . ® Hn ■ Throughout 
the paper, the respective Hilbert spaces are assumed to have 
finite dimensions, Hj = <D dj , j = 1, N. 

One way of solving this problem is to choose a specific ba- 
sis for the Hilbert space and to explicitly parametrize the state 
vectors. This yields a complex polynomial in these param- 
eters, in general of very high order. This is obviously not a 
convex problem in these variables: a solution can be found, al- 
beit not in an efficient manner. For small systems, algorithms 
such as simulated annealing may be employed, delivering up- 



per bounds to the optimal solution, as no control is possible as 
to what extent one is far away from the global optimum. 

The general strategy of this paper is in instances of the 
above type to introduce additional variables, giving rise to one 
vector x 6 R*, x = (xi, Xt) T , which is the objective vari- 
able, parametrizing the product states. The problem is then 
cast into the form of a linear objective function, simply as 

minimize c T x (3) 

with a (fixed) c €E R*, subject to constraints which are polyno- 
mials in the objective variables. These constraints will then be 
relaxed to semi-definite problems. So two types of contraints 
will be encountered in the present paper: 

• Semi-definite constraints: These are constraints of the 
form 

t 

F + 5> S F S >0 (4) 

s=l 

where Fq, F t are Hermitean matrices of arbitrary di- 
mensions. The resulting matrix has to be positive semi- 
definite, therefore, it is referred to as semi-definite con- 
straint. Optimization problems of this type, exhibit- 
ing a linear objective function and semi-definite con- 
straints, are called semi-definite programs fl6ll . Such 
instances of convex optimization problems can be effi- 
ciently solved, for example by means of interior-point 
methods fi^l . Moreover, the idea of Lagrange dual- 
ity JT7I1 readily delivers lower bounds for the problem. 
Typically, the dual optimization problem yields an opti- 
mal value which is identical to the optimal value for the 
primal problem (unless there is a duality gap). Many 
problems in quantum information theory have already 
the form of a semi-definite program IH In fact, it 
may be convincingly argued that to specify the solution 
of a problem in form of a semi-definite program has the 
same status as stating a result in terms of the spectrum 
of a matrix, as this again merely means that efficient 
methods are available to find the eigenvalues of a given 
matrix. 

• Polynomial constraints: This means that we can write 
the constraints as 

gi(x)<o, (5) 

l = 1, L, where gi : R* — > R are real polynomials 
of some degree. Quadratic constraints are of the form 

x T Aix + bjx + a < 0, (6) 

I = 1, L. The matrices A\ are, however, not neces- 
sarily positive semi-definite. This is by no means a mi- 
nor detail: if all matrices A\,..,,Al were positive ma- 
trices, Ai > 0, this would yield a convex quadratic pro- 
gram, which can be efficiently solved (they are in fact 
also instances of semi-definite programs and of second- 
order cone programs). In stark contrast, if the matri- 
ces are not all positive semi-definite, one obtains a very 
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hard, non-convex optimization problem. This structure 
is yet dictated by the problems from quantum informa- 
tion theory at hand. 

The central point is to employ known methods from the the- 
ory of relaxations of non-convex optimization problems, to 
obtain complete hierarchies of cheaply computable approxi- 
mations, approximating the solution as closely as desired. The 
idea of a relaxation is to introduce new variables and to formu- 
late the problem as a convex problem in a larger space. This 
idea can be exemplified in the simplest form of a relaxation, 
the Shor relaxation lUlll . For example, let A\ in Ineq. ^ be a 
matrix which is not positive semi-definite, and let us assume 
that b\ = and c\ = for simplicity. Then, one can still write 
the constraint equivalently as 

tr[XAx] < 0, X = xx T : (7) 

using at x t symmetric matrix X. The equality X = xx T is 
equivalent with the convex constraint 

X > xx T : (8) 

together with the non-convex one X < xx T . Shor's re- 
laxation amounts to taking only the convex part into ac- 
count, thereby delivering an efficiently solvable convex prob- 
lem which yields a lower bound to the original problem full . 
Such relaxations in terms of semi-definite constraints will be 
employed, yet instead of one many such relaxations forming 
a complete hierarchy. 

As pointed out before, we will show that the encountered 
optimization problems can be written as polynomially con- 
strained problems of degree three. That this is possible is 
based on the observation that any Hermitean m x m-matrix O 
for which 

tr[<3 2 ] = 1, tr[0 3 ] = 1 (9) 

is one that satisfies 

tr[0] = 1, O = O 2 , O > 0, (10) 

i.e., it corresponds to a pure state (see also Ref. 0)- This 
follows from the fact that, denoting the decreasingly ordered 
list of eigenvalues of O with (O), the only vector consistent 
with 

m m 

J2^(Of = l, ]TAj(0) 3 = l (11) 

i=l i=l 

is the vector X l (0) = (1,0,. ..,0). The quantities \\{0) 2 
and \\ (O) 3 are unitarily invariant, and hence, the above state- 
ments can be shown to be valid on the level of probability dis- 
tributions. Essentially, X)"=i ^i(0) 2 — 1 already requires all 
absolute values of eigenvalues to be smaller than or equal to 
1, such that the only ordered vector of real numbers consistent 

withEHi Aj(O) 3 = 1 becomes (1,0,. ..,0). 

For systems where the individual constituents are qubit sys- 
tems, dj = 2 for all j = 1,...,N, the constraints can further 
be simplified by merely requiring as constraints tr[Oj = 1, 



tr[0 2 ] = 1, as for Hermitean 2 x 2-matrices these conditions 
alone imply that 

O > 0, = 2 . (12) 

When applied to our specific problems at hand, these con- 
straints will appear in the following form. We will require that 
Hermitean matrices P are, except from normalization, prod- 
ucts of pure states with respect to all constituents. This will 
be incorporated as follows: Denoting with I = {1, TV} the 
index set labeling the subsystems and with tr A:/ - the partial 
trace with respect to all systems except the one with label j, 
the lines 

tr[tr Aj [P] 2 ] = (tr[P]) 2 , (13) 
tr[tr Ai [P] 3 ] = (tr[P]) 3 (14) 

for all j E I indeed enforce that the matrices are products. 
If reductions are pure, the global state must be a pure prod- 
uct state. This can be seen as follows. For states p, the only 
possibility for 

tr[tr Aj [p] 2 ] = l, tr[tr Aj .[p] 3 ] = 1 (15) 

to hold for all j € I is that p is of the form of a product pure 
state, 

p=\cf) 1 )(<j) 1 \® ...®\<j> N )(<j> N \. (16) 

If an additional constant a > is included, these conditions 
read tr[tr Aj - [ap] 2 ] = (tr[ap]) 2 = a 2 and tr[tr/\ 7 -[ap] 3 ] = 
(tr[ap]) 3 = a 3 , which explains the above constraint 11911 . 
Having stated the general strategy, let us now look at the spe- 
cific instances of problems in quantum information we will be 
considering in this paper. 



A. Tests for bi-partite and multi-partite entanglement 

The approach is here to consider for a given state p G 
S(Hi (8> ... <§) TLn) the minimal Hilbert-Schmidt norm with 
respect to the set of separable states. For simplicity of nota- 
tion, we explicitly formulate the optimization problem for the 
instance of full separability, without loss of generality. That 
is, we test whether p can be written as 

n 

»=1 

with {pi]i forming a probability distribution. The question 
whether a state is fully separable is hence equivalent to ask- 
ing whether a state is an element of the convex hull of prod- 
uct vectors with respect to all subsystems. According to 
Caratheodory's theorem l20ll . for any ^-dimensional subset 
S C R m , any point of the convex hull of S can be written as 
a convex combination of at most k + 1 points from S. Hence, 
the number of elements in the convex combination given by 
Eq. d!7l > can be restricted to n = Yl^=i d 2 , again without loss 
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of generality. To decide whether a state p is fully separable or 
not, we may solve the following optimization problem, 

minimize \\p - P\\ 2 2 = tr(p - P) 2 , (18) 
subject to P is fully separable. 

We make use of the Hilbert-Schmidt norm as it is quadratic in 
the matrix entries. 

The task is to write this problem in terms of a polynomially 
constrained problem. Each relaxation (see Section ITTTfc . la- 
beled with h = hmm, h mm + 1, then delivers a lower bound 
of the Hilbert Schmidt distance to the set of fully separable 
states. Hence, asserting that the state is not fully separable 
whenever we obtain a value larger than the one that we accept 
as accuracy of the computation [21], each step delivers a suffi- 
cient criterion for multi-partite entanglement in its own right, 
and the hierarchy is complete in the sense that each entangled 
state is detected by some step. The associated optimization 
problem can now be written as 

minimize x, (19) 
subject to x > tr(p — P) 2 , 

n 

p-J2 p(i) = °> 
i=i 

tr[tr Aj [P«] 2 ] = (tr[pW]) 2 

for all i = 1, n, j G I, 

n 

5>[P»]=1, 

i=l 

P (i) > 0, for alii = l,...,n. 

The line Yl7=i trf-P*- 1 '] = 1 takes the normalization of the 
whole state into account. This is a quadratic program, com- 
bined with a semi-definite constraint for the positivity of the 
matrices PW . As pointed out above, the problem can also be 
formulated as a polynomial problem without a semi-definite 
constraint, but now with constraints that are of degree three. 

minimize x, (20) 
subject to x > tr(p — P) 2 , 

n 

P-Yj P(l) = . 

i=l 

tr[tr Aj [P«] 2 ] = (tr[P«]) 2 , 

for all i = 1, n, j G I. 

tr[tr Ai [P«] 3 ] = (tr[pW]) 3 , 

for all i = 1, ri, j G I. 

n 

5>[P«] = 1. 

i=l 

This is a global optimization problem with polynomial con- 
straints of degree three, but no semi-definite constraint. 

Here one tests the hypothesis that the state is fully separa- 
ble against the alternative that the state is entangled in some 
sense. To assert that the state is multi-particle entangled and 



not separable with respect to any separability class, several 
tests are hence required. In this way, the various classes of 
genuine multi-particle entanglement can be detected. Note 
that when even applied to the bi-partite case, the resulting hi- 
erarchy of semi-definite relaxations is inequivalent to the one 
in Ref. |8|. and also inequivalent to the robust semi-definite 
programming approach in Refs. 0]. The above formulation 
in the optimization problem in terms of full separability still 
does not constitute a restriction of generality, as this includes 
all separability classes with respect to all possible splits. 

Alternatively to the above approach, one may write each 
test in the form of a feasibility problem, a problem with a 
vanishing objective function, 

minimize 0, (21) 
subject to p satisfies the test of step 

h = h m { n , ft-min+i, ••• in the hierarchy. 

Either one finds no solution (which is to say, the problem is 
not primal feasible), and one can assert that the state is not 
fully separable, or one has to go on one step in the hierarchy. 
In each step of the hierarchy, forming a semi-definite problem, 
the dual problem can then be employed to prove the infeasi- 
bility of the above primal problem serving as a certificate fl^l 
(see also Ref. @]). 

In general the total problem can in each step be written as a 
semi-definite problem of the form 

minimize 0, (22) 

T 

subject to Hq + z s H s > 0, 

s=l 

with appropriate matrices H s , s = l,...,T. The associated 
Lagrange dual problem fl7ll is again a semi-definite program, 

maximize —tr[ZHo] , (23) 
subject to tr[ZH s ] =0,s= 1, T, 
Z>0. 

In the context of our feasibility problem above, any feasible 
solution of the dual problem with tr[ZPo] < proves the in- 
feasibility of the primal (original) problem. That is, we can 
use the dual problem to prove properties of our original prob- 
lem at hand. 

Finally, it is important to point out that in problem dl9> . one 
may keep the semi-definite constraint provided by the last line, 
and look for the intersection of the feasible sets of the semi- 
definite part and the constraint set of the relaxations. Then, 
in each step we can either assert whether the state is entan- 
gled, or one cannot say whether it is entangled or not. In this 
way, it may happen, yet, that the state is entangled, although 
this entanglement is not detected in any step of the hierarchy. 
One hence obtains a hierarchy of sufficient criteria, albeit one 
which is not necessarily asymptotically complete. This will 
be discussed in more detail in the section on the hierarchy of 
relaxations. 

In an implementation of this optimization problem, one 
has to choose a basis of Hermitean matrices for each Hilbert 
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space, 

<rq}, (24) 

for j = 1, N, suppressing an additional index labelling the 
subsystems. These Hermitean matrices satisfy tr[<7i] = 1 and 

tr[a fc ]=0, k = 2,...,d 2 j , (25) 

and have a Hilbert-Schmidt scalar product 

tr[cr fc er/] = S.djSki (26) 

with a dimension dependent constant (and similarly for 
terms of third order). For the case of qubit subsystems, the 
appropriately normalized familiar Pauli matrices can be taken. 
In terms of this basis of Hermitean matrices, the matrices pW 
and P can be written as 

p {i) = E < 2? ) 
p = E (28) 

K=(ki,...,k N ) 

where k = (fci, fc^r), is a multi-index, with fcj = 1, ef| 
for j e /, and 

S K = cr fcl ® crfe 2 ® ... <g> a kN . (29) 

This parametrization will be used in the section presenting nu- 
merical examples. Before we present the hierarchy of relax- 
ations explicitly, we discuss the other applications which are 
similar in structure from the point of view taken in this paper. 

B. Non-decomposable witnesses 

Optimization problems of the type of the one in Eq. 10 of- 
ten appear in the construction of entanglement witnesses 12311 . 
An entanglement witness is a Hermitean observable W = W' 
with the property that tr[Wp] > holds for all separable p, 
thus a negative expectation value signals the presence of en- 
tanglement. So entanglement witnesses can be used for an 
experimental verification that a given state is entangled, and, 
in fact, they have already been implemented 12411 . 

The detection of entanglement is not only of interest for 
fundamental reasons, it can also be of practical interest. This 
is the case in quantum crytography, since it has been shown 
that the provable presence of quantum correlations in such 
protocols is a necessary precondition for secure key distilla- 
tion l25ll . Entanglement witnesses are particularly suited to 
deliver this entanglement proof, even when the quantum state 
shared by the users cannot be completely reconstructed. In 
turn, by measuring all accessible witnesses, one can decide 
whether the measurable correlations of the state origin from an 
entangled state or may be compatible with a separable state. 

There are many strategies to construct entanglement wit- 
nesses I2H l27ll . As an example in which such optimi- 
sation problems occur we choose the construction of non- 
decomposable witnesses for PPT entangled states |27]. These 



are entangled states which have a positive partial transpose 
I2H1 . We discuss our example in the bi-partite setting for 
simplicity. In the theory of PPT entangled states the extreme 
points of the set are of central interest, then often referred to 
as edge states. A state p is a PPT entangled edge state if it has 
a positive partial transpose, while for all product vectors \a, b) 
in the range of p the vector |q, b*) is not in the range of the 
partially transposed p Ta I23l 12711 . Here * refers to complex 
conjugation. 

To construct witnesses for these states, one proceeds as fol- 
lows: Let R = K{p) and Q = K(p TB ) be the projectors 
onto the kernels of p and of p TB , respectively. Then a witness 
allowing the detection of the state p is given by l27ll 

W' = R + Q Tb - el (30) 

where 

e = mintr[|a,6)(a,6|(P + Q Tf3 )]. (31) 

\a,b) 

Since p is an edge state, we have e > 0. This implies that 
tr[M^'p] < 0, thus, p is detected. Also, it is clear that the diffi- 
cult part of this construction is the minimization procedure in 
Eq. ( 13 11 1 - which is just of the type of Eq. 0. 

This method can also be used to obtain a finer witness from 
a given one W, i.e., a witness that detects the same states as W 
- and more. If e = minuM tr[|a, b) (a, b\W] > then W — el 
is a finer witness than W. This can also be applied in the sce- 
nario where only a restricted set of observables is available, 
since the observable 1 is always accessible. In practical sit- 
uations, given a particular implementation of a quantum key 
distribution (QKD) scheme, it is sufficient to obtain one rele- 
vant entanglement witness as a first step towards the demon- 
stration of the feasibility of the scheme. Here is where the 
method presented in this section can be used. Although this 
method requires, as a starting point, to have already a valid 
entanglement witness for the given QKD protocol, note that 
this operator does not need to be an entanglement witness in 
the strict sense, but can be a positive operator from the re- 
stricted set which is more easy to characterize than an entan- 
glement witnesses l3lll . Moreover, during several steps of the 
method, better entanglement witnesses can be obtained from 
it, belonging to the same restricted set. 

The construction of the witness above can also be used for 
multi -partite PPT witnesses l32il . The other approaches for 
the construction of entanglement witnesses also need simi- 
lar optimization processes 12611 . For the sake of generality, 
we formulate the optimization strategies directly in the multi- 
partite setting as in Eq. (|2ji. For a given entanglement witness 
W = W\ the optimization problem looks as follows: the aim 
is to solve the problem 

minimize x, (32) 

subject to x>tr[WP], 

tr[tr Ai [P] 2 ] = 1, forallje/, 
tr[tr Aj [P] 3 ] = 1, for all j e I. 

Again, for multi-party qubit systems this can be written as 
a polynomially constrained problem with polynomials of de- 
gree two. 
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C. Estimating the geometric entanglement to quantify 
multi-particle entanglement 

The same tools can be used in order to quantify multi- 
particle entanglement for pure quantum states. Needless to 
say, the question of quantifying multi-particle entanglement 
is much more involved that the analogous question in the bi- 
partite setting: in the bi-partite setting, the degree of entan- 
glement of pure states can be uniquely quantified in terms of 
the entropy of entanglement. Any pure state can be asymp- 
totically reversibly transformed into any other, the achievable 
rates being given by just this measure of entanglement. In this 
sense, any bi-partite entanglement of pure states is essentially 
equivalent to that of the maximally entangled pair of qubits, 
which forms the so-called minimal reversible entanglement- 
generating set (MREGS) l33ll . The situation is very differ- 
ent in the multi-partite case, where the MREGS have not even 
been identified for three-qubit systems, let alone for more gen- 
eral settings l34ll . In the view of this fact, several more prag- 
matic (and inequivalent) measures of entanglement have been 
proposed, reasonably grasping the degree of multi-particle en- 
tanglement 1 3 5113a. 13711. To evaluate these quantities typically 
amounts to solving a computationally hard problem. 

One of the reasonable quantities to quantify multi-particle 
entanglement is the geometric measure of entanglement l35l 
I36H38II : for a given state vector \tp) <E Hi <8> ... (£> TLn, essen- 
tially, entanglement is then quantified in terms of the solution 
of the maximization problem 

A 2 = max|(?/>|0)| 2 , (33) 

such that the geometric measure of entanglement becomes 

£(|V)(V'I) = 1-A 2 . (34) 

The maximization is performed over all state vectors \<p) 
which are products with respect to all subsystems. Setting 

p = \ip) (ip\ and 

p = 10X01 = \<f>i)(4>i\ ® ... ® \4>n)(4>n\, 05) 

we arrive at 

minimize t, (36) 

subject to tr[Pp] + t > 1, 

tr[tr Aj [P] 2 ] = 1, for all j g I, 
tr[tr A .,[P] 3 ] = 1, for all j g I, 

which is the same optimization as in the previous subsection, 
except from one line in the list of constraints. 

D. Entanglement witnesses based on second moments and 
entropic criteria 

In this subsection we will consider again entanglement wit- 
nesses, but not in the original sense, which involve only ex- 
pectation values of Hermitean operators. It is also possible to 
introduce nonlinear functionals with similar properties: these 



are entanglement witnesses based on second moments, on 
variances of observables. Such entanglement criteria based 
on second moments are very popular in the study of infinite- 
dimensional quantum systems having canonical coordinates 
1I39IL There, to measure arbitrary observables is often by 
far unfeasible, whereas the estimation of second moments of 
canonical coordinates is very accessible. In optical systems, 
the appropriate measurements are available in homodyne de- 
tection. Similarly, one may also for finite-dimensional sys- 
tems look at variances rather than at first moments themselves 
lEoll . In Ref. BUl . second order witnesses related to variances 
of operators were constructed, and the relation to entangle- 
ment criteria for continuous variable systems based on second 
moments was shown. The advantage in a practical context is 
that one specifies some observables which are the most ac- 
cessible, and tests whether the obtained second moments are 
consistent with a separable state. The application of such a 
test always requires the solution of an optimization problem 
as pointed out below. 

Let us specify a set of observables Mi, Mr-. Then we 
can define the real symmetric K x K covariance matrix 7 P of 
a state p associated with these observables as 



(7pk/ = HM k Mip]+\x[MiM k p})/2 

- tr[M fc p]tr[Af,p], (37) 



with k, I = 1,...,K. This is completely analogous to the 
familiar covariance matrix of systems with canonical coordi- 
nates. Then, it turns out that - in the previous notation - any 
fully separable state p has the property that there exist states 



Pi , —,Pni i = -> n ' ( 38 ) 



n = Yij=i dj, and probability distributions {pi}i such that 



-Yp>Y,MpP9...9p$' (39) 



So one would fix those observables that are the most acces- 
sible, and estimate the appropriate second moments. This 
would yield an estimate of the elements of the covariance 
matrix to some accuracy. Then, the question that arises is: 
do states and probability distributions exist that satisfy Ineq. 
j39t? If not, we can conclude that the state must have been 
entangled. It is important to note that this judgement is not 
based on the knowledge of the entire state, but only on the 
knowledge of the covariance matrix with respect to a previ- 
ously selected set of observables. This is a problem that can 
be cast into a feasibility problem, again in the form that we 
envision. As it is a feasibility problem, the objective function 
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can be set to zero. This can be written as follows, 

minimize (40) 

n 

subject to Q - = 0, 

i=l 

Qg = (tr[M fe MiP«] + tr[M;M fe P«])/2 

-tr[M fe P«]tr[AfjPW], 
tr[tr Aj [p(''] 2 ] = (tr[P«]) 2 

for all i = 1, n, j e 7, 
tr[tr Aj [pW] 3 ] = (tr[p( i )]) 3 

for all % = 1, n, j G /, 

n 

5>[P«] = 1, 

7 P - Q > o. 

In each step of the hierarchies of relaxations, we can assess 
whether there is a feasible solution or not. If there is no fea- 
sible solution in some step, we can conclude that the state is 
entangled, and multi-particle entanglement is hence detected. 
This is now a problem which is still a combination of a poly- 
nomially constrained problem, together with a semi-definite 
constraint. Here, two strategies may be applied: 

On the one hand, one can keep the semi-definite constraint, 
and can proceed as pointed out in Subsection A. This means 
that one can in each step assert that the state was entangled, 
or one has to go one step further. This is computationally 
cheaper, but comes at the price of losing asymptotic complete- 
ness. For practical purposes, however, this method is expected 
to be the method of choice; in particular in the light of the 
fact that for a given set of observables, typically not every en- 
tangled state is anyway detected by the entanglement witness 
based on second moments. 

The other strategy, on the other hand, is to formulate 
7p — Q > as a set of polynomial constraints. A Hermitean 
matrix is positive if and only if the determinants of all its sub- 
matrices are positive (see also Ref. 13111 1 This gives rise to 
a set of polynomial constraints, for which the relaxations can 
be applied, leading to an asymptotically complete hierarchy of 
tests. This comes at the price of being computationally more 
expensive. 

In view of Refs. l40H4lll . it is useful to employ arguments 
along the following line: if we can show that no fully sepa- 
rable state can have the image that we estimate in an experi- 
ment, we can assert that the state must have been entangled. 
This observation can, while being fairly obvious, still be prac- 
tically very relevant. For example, we may for any observable 
M = Aft look at the minimum of second moments that are 
consistent with a separable state, i.e., the solution of 

minimize tr[M 2 P] - tr[A/P] 2 , (41) 
subject to P is fully separable, (42) 



gives rise to the optimization problem, 

minimize x, (43) 
subject to x > tr[M 2 P] - tr[A/P] 2 , 

n 

p-Yj p(l) = °> 

i=l 

tr[tr Aj [P«] 2 ] = (tr[P«]) 2 , 

for all i = 1, n, j € I. 

tr[tr Aj [P«] 3 ] = (tr[P«]) 3 , 

for all i = 1, n, j € I. 

f>[P«] = l. 

i=l 

It should be clear at this point that the same method can be 
used for entanglement criteria based on linear entropies, that 
is, p-norms forp = 2 (see, in the rich literature on the subject, 
e.g., Refs. l42ll ). For any expression that is linear in the linear 
entropies of the whole state p, 

l|p||2=tr[ (0 2 ] (44) 

of a multi-partite system and in the linear entropies of the re- 
ductions 

lltr^/Wlb =tr[(tr jV [p]) 2 ], j = l,...,n, (45) 

one can in the same manner find the largest value consistent 
with a separable state. Any state that delivers a larger value is 
then clearly entangled. In practical considerations, these lin- 
ear entropies can be estimated in a fairly feasible manner l42ll . 
for example when assessing entanglement in Bose-Hubbard- 
type models. 

E. Maximal output purities of quantum channels 

Similar arguments, it shall finally be briefly discussed, can 
immediately be applied to assess minimal output purities of 
channels 

k 

p^e(p) = Y,RipRl (46) 

i=l 

where X)i=i ^JPi = 1, with respect to p-norms for p = 2 
(and other integer p), see, e.g., Refs. ll43ll . One may then in- 
vestigate the maximal output purity 

is 2 (£) = max \\S(p)\\ 2 , (47) 
p 

where it does not constitute a restriction of generality to maxi- 
mize not over all states p, but merely over all pure states (com- 
pare also Ref. 0]). The central question here is to see whether 
this quantity is multiplicative in general. That is whether gen- 
erally 



and use this as a criterion for detecting entangled states. This 



V2(£i®e 2 ) = V2(£ 1 )u2(£ 2 ) (48) 
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holds, which means that it is never an advantage to allow for 
entangled inputs when maximizing the output purity. In the 
previously used language, this optimization problem can be 
written as follows, 

maximize x, (49) 

k k 

subject to x < tr[(^ i? l Pi?])(^ RjPR.])], 

tr[P 2 ] = 1, 
tr[P 3 ] = 1, 

as a polynomially constrained problem with polynomials of 
degree three. 

III. COMPLETE HIERARCHIES OF RELAXATIONS TO 
APPROXIMATE THE SOLUTIONS 

We will now state how the theory of relaxations can be ap- 
plied to the described problems in entanglement theory. In all 
of the above cases (except in Subsection D), we obtained an 
optimization problem of the following structure: for x G R*, 

minimize c T x (50) 
subject to gi(x) > 0, 

for / = 1, L, the global optimum value being denoted as 
p* , where gi is a polynomial of at most degree three. The 
constraint set given by 

M = {x e R* :gi(x) > 0, I = 1,...,L} (51) 

is not a convex set. We may however apply Lasserre's method 
of semi-definite relaxations to treat this part, see Appendix 
A. This will yield a sequence of semi-definite programs, la- 
beled with an index h = h m i n , h m i n + 1, such that each 
of the efficiently solvable steps yields an approximation of the 
original problem. The minimal step h mm is 1 if the highest 
degree of the constraint polynomials is 2 and h m m = 2 if con- 
straints of degree 3 are required. The case h = /i m ; n is the 
first semi-definite relaxation leading to the first approxima- 
tion, h = hmin + 1 is the second, and so on. Often, in practice 
the global optimum is already achieved after a small number 
of steps in the hierarchy. 

A. Semi-definite relaxations 

Instead of considering an optimization problem in x G R*, 
this is turned into an optimization problem in a larger real vec- 
tor y G R 2h , D^h being a natural number defined in Ap- 
pendix A. This larger dimension is due to the uplifting proce- 
dure used in Lasserre's method for approximating the 
quadratic part of our problems, and goes back to work in Ref. 
14411 . For each instance of the hierarchy of semi-definite pro- 
grams, the objective function will be the same, but uplifted, 
namely 

V — > d T y, (52) 



where 

d T = (0,ci,...,c t ,0,...,0), (53) 

with c G R* being defined as above. Lasserre's method now 
gives rise to a sequence of semi-definite programs approxi- 
mating the solutions of 

minimize c T x (54) 
subject to gi{x) > 0, 

for I = 1 , . . . , L in the following form: For h = h m - m , h m [ n + 
1, each instance is of the form 

minimize d T y, (55) 
subject to F [h] (y) > 0, 

G[ h] (y)>0, l = l,...,L 

with matrices F^ h \y) and G{ h '(y) that are linear in the ele- 
ments of y that increase in dimension with increasing h (see 
Appendix A). This method is based on recent results in real 
algebraic geometry, see also Ref. fl5ll . 

In Ref. lll3ll convergence to the solution of j54t is guar- 
anteed if certain conditions are satisfied: Convergence in 
the limit h — > 00 is guaranteed if there exist polynomials, 
u , ui, ul, all sums of squares, such that the set 

L 

{x G R* : uo(x) +^^ui(x)gi(x) > 0} (56) 

is compact. This is, however, the case in all of the specific 
situations from entanglement theory considered above. The 
set in Eq. d56t is compact if there exists an I € {1, L} such 
that the set 

{x G R* : gi(x) > 0} (57) 

is compact. In each of the discussed cases, we find that due 
to the linear constraints incorporating the trace requirement 
and the quadratic constraints coming from the purity of the 
reduced states, there exists an a > such that a 2 — \\x\\ 2 > 
for all x G M.. This follows from the fact that for each of the 
involved matrices, the trace is bounded from above, and pos- 
itivity of the matrices enforces boundedness of all elements. 
Hence, to ensure asymptotic completeness, we may add the 
constraint g^+itx) = a 2 — ||x|| 2 > to the list of quadratic 
constraints, such that the condition in Eq. ( 15 1\ is certainly sat- 
isfied. Hence, one can conclude that 

min d T y — ► min c T x (58) 

for h — > 00, and for 

M [h] = {ye M, D2h : F^(y) > 0, 

G\ h] (y)>0, l = l,...,L + l}. (59) 

This is not only meant as a numerical procedure: instead, as 
each step is an analytically accessible semi-definite program, 
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in each step one may assess the approximations with analyti- 
cal means. Moreover, symmetries of the involved states under 
certain groups can be carried over to symmetries in the Her- 
mitean matrices in the semi-definite programs, similarly to the 
strategy employed in Ref. |5] for semi-definite programs, and 
in Ref. fl| for convex but not semi-definite programs. 

In Subsection D of Section|II] we encountered an additional 
semi-definite constraint. Then, Lasserre's method may be ap- 
plied using polynomials of higher order, as described above. 
Or, one may combine the semi-definite relaxations with the 
semi-definite constraint itself. This gives rise to a hierarchy of 
sufficient tests, without the property of asymptotic complete- 
ness. To see how they can be combined, let us consider an 
additional semi-definite constraint such as j p — Q > 0. In 
terms of the y G R D2h , we have the feasible set of the addi- 
tional semi-definite constraint 



For example, 



T = {y £ M, D2h : F 



Vs+iF s > 0}, 



(60) 



with appropriate matrices Fq, F t . Therefore, we can write 
the full hierarchy of semi-definite programs as 

minimize d T y, (61) 
subject to F [h] (y) > 0, 

G|%)>0, l = l,...,L 
t 

F + J2ys+iF s >0, 

s=l 

h = h m [ n , h m i n+ i, ... being the label of the element of the hier- 
archy. The projection of the feasible sets M\ h ^ onto the plane 
of first order moments, i.e., onto the plane 



{neK° !k :y=(0,y 2 ,...,yt+u0,-,0)}, 



(62) 



conceived as a subset of IR', converges (pointwise) to the con- 
vex hull of M Therefore, we have that 



mm a y < p 



(63) 



for all h — > oo. Moreover, iuiiiy^zj^[h]f^jr d y is a monotone 
increasing sequence in h, such that the sufficient criteria be- 
come more powerful with an increasing order of the hierarchy. 



B. Size of the relaxations 

A relevant issue is how large the semi-definite relaxations 
are in each step of the hierarchy. In the worst-case scenario, 
where the polynomial constraints is a polynomial involving 
all basis elements of the basis of polynomials of the respec- 
tive degree, one obtains the subsequent sizes of the relaxation 
matrices. The matrix F<- h ' is of dimension Dy L x Dy t (for a 
definition of Dh, see Appendix A). As a formula for Dh we 
arrive at 



E 

fc=0 



t + k-1 

k 



(64) 



D 2 = 1 + t + 



t(t+l) 



(65) 



In the number of variables t, this is a manifestly polynomial 
expression. In step h the vector y is of the length D^h- No- 
tably, in each of the steps, the effort of a numerical solution 
of the associated semi-definite program is polynomial in the 
dimension of the matrices fl^l . Hence, each problem can be 
solved in an efficient manner. 

In terms of the step h in the hierarchy, it turns out that the 
scaling is also polynomial. Approximating the above sum by 
an integral expression, we arrive at 



0{h'). 



(66) 



That is, for a fixed number of variables (which is the setting 
considered here), the size of the vector of the objective vari- 
ables increases also only polynomially in the step h in the 
hierarchy. Moreover, in many small and medium size prob- 
lems, the programm detects optimal solutions in the first iter- 
ation steps at relatively low computational cost 14511 . Also, the 
sparsity of the moment matrices may be exploited. The issue 
of computational effort will be discussed in more detail else- 
where. Another point of interest is that it is possible in some 
cases to trade in a lower number of variables t for a higher 
lowest relaxation step /i m ; n , as in the examples in the subse- 
quent section. In some cases, this might simplify the problem, 
as in our example in the next section. 



IV. NUMERICAL EXAMPLES 

In this section we present some numerical examples, in or- 
der to show that the approach is also feasible in practice. We 
will provide three examples, two for the geometric measure 
of entanglement and one for the construction of entanglement 
witnesses for bound entangled three-qubit states. 



A. Geometric measure for three-qubit states 

Let us start with the calculation of the geometric measure 
of entanglement for three-qubit states. As we have shown in 
Section||l]the computation of the geometric measure of entan- 
glement for a given pure three-qubit state vector \tp) requires 
essentially the calculation of 

A 2 = max | (a, b, c\ip}\ 2 

\a,b,c) 



(67) 



As already mentioned above, we use here a different 
parametrization than the general one described in Section II. 
In terms of the Pauli matrices forming a basis of Hermitean 
matrices, we can write 

1 3 

I^K^I = g E ^ijk(<ri®o-j®a k ), (68) 

i,j,k=0 
1 3 

\a,b,c)(a,b,c\ = - ^ cubjCkfa <g> (7j ® cr k ), (69) 



i,j,k=0 
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where Aooo = a o = &o = c o = 1- The coefficients Xijk, 
k = 0, 3, are determined from the known state vec- 
tor We have to impose constraints that guarantee that 
Pa is a pure state on the coefficients (01,02,03) describing 

the state pa = Y^i=o a i a i/^ C 311 ^ similarly (&i,&2>&3) an d 
(ci, C2, C3)). We have seen before that for qubit systems, in- 
stead of requiring tr[p A ] = 1 and tr[p A ] = 1, we may alterna- 
tively merely require that 1i[pa] = 1 an d trfp^] = 1 (where 
tr[jOA] = 1 is already a consequence of the parametrization). 
So we arrive at the optimization problem 



minimize,. 



1 3 



Cfc, 



(70) 



i,j^k— 



subject to 



+ &1 + 6| = 1, 
c? + c| + c§ = l. 



This polynomial optimization problem can be solved with the 
help of Lasserre's method, see Appendix A. For these calcu- 
lations the package GloptiPoly |Ej based on SeDuMi Eclll 
is freely available, and we have used it for our calculations. 
The package GloptiPoly has a number of desirable features, 
in particular, it provides a certificate for global optimality. 

Note that with this parametrization, the number of objective 
variables is AN, N being the number of qubits, in contrast to 
A N parameters which are necessary to parametrize a general 
N qubit state as decribed in Section II. A. From Eq. i70\ it is 
clear that the objective function will be a polynomial of degree 
N which increases h m i n , see Appendix A. 

First, we present a nontrivial example for the calculation of 
the geometric measure of entanglement, in a case where its 
value is already known. In this way we can test our methods. 
We aim at computing the geometric measure of entanglement 
for state vectors of the form 



(71) 



s G [0, 1], where \W) and \W) are state vectors of three-qubit 
W states in different bases, 



\W) 
\W) 



(|001) 
(|011) 



|010) 
|101) 



|100»/V3, 
|110))/V3. 



(72) 
(73) 



For the geometric measure of entanglement of \ip(s)) a for- 
mula has been developed in Ref. ll36ll . exploiting the permuta- 
tion symmetry of the states. The comparison between the the- 
oretical value and the numerical calculation using Lasserre's 
method for ft, = 2 is shown in Fig. [2 Details of the perfor- 
mance are summarized in Table|I] The results indicate clearly 
the usefulness of the presented approach. As a matter of fact, 
this is a case where already a very small number of steps in 
the hierarchy detects the global optimum, as is typical for this 
relaxation, as has been pointed out in Ref. ll45ll . based on nu- 
merical experiments. 




FIG. 1: The numerical values of the geometric measure of entan- 
glement E of the family of states of Eq. 17 1 i . plotted on top of the 
analytical values of Ref. l3(Hl . 



Subsection 


Relaxation h 


# variables 


clim(y) 


CPU time 


A 


2 


9 


714 


10.92 s 


B 


2 


12 


1819 


103.97 s 


C 


2 


9 


714 


6.14 s 



TABLE I: Details of the relaxations in the three numerical examples 
discussed above for one point of each example. The provided CPU 
time refers to a machine with a Intel Xeon Processor, 2.2 GHz, 1GB 
Ram, using GloptiPoly 2.2e 0], SeDuMi 1.05 and MatLab 
6.5.1.199709 (release 13). In all cases h = /i m in = 2, so that the 
result was obtained after the first relaxation step. 



B. Geometric measure for 4-qubit states 

We calculate the geometric measure of entanglement also 
for the following one parameter family of state vectors 



\Mp)) = v^lGHZ') - v / T^PlV' + ) ® IV> H 



where 



|GHZ') = (10011) + |1100))/V2, 



(74) 



(75) 



|V' + ) = (|01) + |10))/a/2, andp e [0,1]. The state vec- 
tor |^4(2/3)) corresponds to the 4-qubit singlet state, i.e., the 
state vector satisfying 



U 



84 



(76) 



for all unitary U 14^1 . For the two individual states in the 
above superpositions in Eq. il 4> . the geometric measure can 
be directly evaluated 13611 : For p = 1 we find A 2 = 1/2, and 
for p = we obtain A 2 = 1/4 from A 2 + = 1/2. The nu- 
merical results for the geometric measure of entanglement for 
other values of p are plotted in Fig. [2] It is interesting to note 
that at the singlet value p = 2/3, the behavior of the geomet- 
ric measure changes. From there up to p = 1 the optimum is 
attained for the choices |0011) or 1 1 100) of the product state 
which gives rise to the linear behavior. 
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FIG. 2: The numerical values of the geometric measure of entangle- 
ment E of the family of states of Eq. I74i . 



product states has to be performed with respect to |3J 
W = -(|000}(000| + |H1)(111| 



1 



1 + c 2 
1 

T+h 2 
i 



1 + a 2 

1 c 

- -I 

2 



c 2 |100)(100| + |011) <011 1 
|010)(010| + 6 2 |101)(101| 
|001)(001| + a 2 |110)(110 



(|000) 



1 + c 2 
1111 



1 + b 2 l + a s 

|iii)(ooof 



(79) 



The numerical results are plotted in Fig. (|3jl. Again, the global 
optimum is achieved, and the found values agree with the val- 
ues found in Ref. 13211 . For details concerning the relaxations, 
see Table U 



The family of states specified in Eq. i74i is invariant un- 
der the exchange (AB) <-» {CD). Because of this sym- 
metry, one may without loss of generality assume that the 
product state vector leading to the maximal value of A 2 is 



given by 



where <£i 2 ) 



cos0i, 2 |O) 



="71,2 



sin6*i 2|l), where the optimal phases can be shown to 
be xi — X2 — Vi = V2 = 0. This gives rises to an op- 
timization problem with polynomial constraints with only 4 
variables which can be solved exactly by GloptiPoly. The re- 
sults coincide with the results above. 



C. Witness for 3-qubit PPT entangled states 

Employing the same strategy, we would like to calculate the 
value of e as defined in Section HTIfor the family of witnesses 
constructed for the PPT (bound) entangled states 



p = (a|001)(001| + b|010)(010| + c|011)(011| 

+ -|100)(100| + -|101)(101| + -|110)(110| 

c b a 



(77) 



2|GHZ)(GHZ| )/{2 + a + b + c + l/a+ 1/6+ 1/c) 



where a, b, c > 0, ab ^ c, and 



|GHZ) = (| 000) + |111 



(78) 



In Ref. [ 32] , upper bounds for the values of e were obtained by 
using a multi-variable minimization routine for the parameter 
range a = b = 1/c € (0,1). The minimization over the 
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FIG. 3: The numerical values of e for W of Eq. 1791 plotted on top 
of the results of Ref. 1321. 



V. SUMMARY AND OUTLOOK 

In this paper, we have considered several problems in en- 
tanglement theory with the tools and language of convex opti- 
mization. The central point of the paper was that many prob- 
lems, where a minimization over pure product vectors is re- 
quired, can be written as instances of certain optimization 
problems involving polynomial constraints of degree two or 
three, or with additional semi-definite constraints. For such 
polynomially constrained problems, which are generally in- 
stances of non-convex optimization problems, hierarchies of 
semi-definite relaxations can be found. In this sense, one addi- 
tional intention of this paper is to communicate these recently 
achieved results in the theory of relaxations and to show that 
they can be fruitfully applied in the quantum information con- 
text. One arrives at hierarchies of more and more refined tests 
detecting entangled or separable states, or better and better 
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lower bounds to optimization problems. In all instances, re- 
cently achieved known results from semi-algebraic geometry 
guarantee that asymptotically, the achieved minimum is in- 
deed approaching the globally optimal one. In this sense, the 
statements are similar in spirit with yet more versatile than 
the ones presented in Refs. 0] . Moreover, we have seen that 
the size of the optimization problems to be solved in each test 
grows polynomially with the steps in the hierarchy, and that 
for small problems, often already a small number of steps is 
required to find the exact solution. 

The presented method is on the one hand meant as a nu- 
merical method to achieve good bounds to problems that are 
of relevance in the study of multi-particle entanglement, in 
the construction of entanglement witnesses in the bi-partite 
and multi-partite case, in the context of quantum key distri- 
bution, and to assess maximal output purities. On the other 
hand, each instance of the hierarchy delivers a semi-definite 
program which is readily accessible with analytical methods, 
and where properties of the Lagrange-dual can be exploited. It 
is hoped that these techniques shed new light on the structure 
of optimization problems underlying the questions of entan- 
glement and separability of several constituents. 

Finally - and shifting perspective to some extent - it seems 
worth noticing that very similar techniques may be expected 
to be useful tools to assess ground state properties of many- 
body Hamiltonians. Often variational approaches deliver al- 
ready a good approximation to properties of the true ground 
state. For example, in the Gutzwiller ansatz for the ground 
state of the Bose-Hubbard model in a lattice one optimizes 
the energy functional over product states with respect to the 
sites. Similar techniques can be used for spin systems and ex- 
pressions for ground states in terms of matrix product states. 
Then, relaxations in the way discussed above could poten- 
tially be applied for a reasonable size of the system. Such 
studies could complement numerical techniques yielding up- 
per bounds, such as simulated annealing techniques, deliver- 
ing provable lower bounds for the ground state energy. 
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VII. APPENDIX A: LASSERRE'S METHOD 

For completeness, in this appendix we briefly sketch the 
method to construct sequences of semi-definite relaxations of 
global optimization problems with multi-variate real-valued 
poly nomial objective function and constraints due to Lasserre 
IU3II . The class of problems is of the following form: 

minimize c T x, (80) 
subject to gi{x) > 0, l = l,...,L, (81) 

where gi, gt/, : IR* — > 1R are real-valued polynomials of 
degree two or three. Although we consider only polynomials 
of degree of at most three, it will be convenient to formulate 
the subsequent sequence of semi-definite programs in terms 
that formally involve higher-order polynomials. For any r £ 
IN, we consider the basis of polynomials of degree r in the 
variables x\,...,Xt as 

{l;x\, ...,x t ;xl,xix 2 , ...,x 1 x t ;x\,x 2 xz, ...,x r t ), (82) 

in this ordering. The dimension of this basis will be denoted 
as D r . For clarity of notation, we will not specify t as an 
index, as this will stay the same throughout the procedure. 
Any polynomial of degree of at most r can then be identified 
with a vector p £ H Dr . It is convenient to introduce two 
labelings, connected with each other by a function 

f r : {1, ...,D r } -> I a = (ai, ...,a t ) :^]a s <r , (83) 
^ s=l ' 

such that the i-th element z, i = 1, D r , of the basis given 
by Eq. d82l > is written as 

/ 

z = l[xf, (84) 

?:=i 

characterized by a — (ai,...,a t ) £ Nq. Note that for a 
given k £ N there are ( t+k k ~ 1 ) possible vectors a such that 
a s = k. It follows that the dimensions Dh are given by 
Eq. (|6l. 

In the following we give the required matrices from 
Lasserre's method for general polynomials Il3ll and discuss 
the cases occuring in the paper explicitly afterwards. Let Si 
be the degree of the polynomial constraint I £ {1.....L} 
and f"<5//2] be the smallest integer greater than or equal to 
Si/2. We assume that the objective function is linear, which is 
no restriction of generality, as other polynomials can always 
be incorporated in the constraints as in Section II. Then the 
first possible relaxation step of Lasserre's method is h m i n = 
max; |"<$z/2] . For h > h m i n the matrix (y) is of dimension 
Dh x Dh and linear in a vector y £ R D2h , 

[ F[ftl ^i,i=^(/4i) + M))' < 85 > 

In turn, the matrices (y), one for each of the constraint 
polynomials, I = 1, L, are of dimension Dj x D^ , where 
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hi = h — \5i/2~\. Each polynomial gi is characterized ac- 
cording to the above procedure by a vector The matrices 

Grj (y) are then defined as 

[cfdOiw = E% 1 (a)%- hi (f- hl V)+f h (j))+<*y (86 > 



For qubits, /i m i n = 1, because the maximal degree of the con- 
straint polynomials is 2. For higher dimensional systems, the 
highest occuring order is 3 due to the positivity constraints. In 
this case, ft, m ; n = 2. 



[1] L. Gurvits, Annual ACM Symposium on Theory of Comput- 
ing, Proceedings of the thirty-fifth ACM symposium on theory 
of computing, San Diego, CA, USA, ACM Press, June 9-11 
(2003). 

[2] R.F. Werner, Phys. Rev. A 40, 4277 (1989). 

[3] S.L. Woronowicz, Rep. Math. Phys. 10, 165 (1976); A. Peres, 
Phys. Rev. Lett. 77, 1413 (1996); M. Horodecki, P. Horodecki, 
and R. Horodecki, Phys. Lett. A 223, 1 (1996); M. Lewenstein, 
D. BruG, J.I. Cirac, B. Kraus, M. Kus, J. Samsonowicz, A. San- 
pera, andR. Tarrach, J. Mod. Opt. 47, 2481 (2000); O. Rudolph, 
J. Phys. A 33, 3951 (2000); M. Horodecki, P. Horodecki, and 
R. Horodecki, in Quantum information: An introduction to ba- 
sic theoretical concepts and experiments, ed. G. Alber et al. 
(Springer, Heidelberg, 2001), pp. 151; B.M. Terhal, Theor. 
Comput. Sci. 287, 313 (2002); K. Eckert, O. Guhne, F. Hulpke, 
P. Hyllus, J. Korbicz, J. Mompart, D. BruG, M. Lewenstein, and 

A. Sanpera, Quantum Information Processing, G. Leuchs and 
T. Beth (Eds) (Wiley- VCH, Weinheim, 2003). 

[4] K. Audenaert, J. Eisert, E. Jane, M.B. Plenio, S. Virmani, and 

B. de Moor, Phys. Rev. Lett. 87, 217902 (2001). 

[5] E.M. Rains, IEEE Trans. Inf. Theory 47, 2921 (2001). 

[6] M. Jezek, J. Rehacek, and J. Fiurasek, Phys. Rev. A 65, 060301 
(2002); K. Audenaert and B. De Moor, Phys. Rev. A 65, 030302 
(2002); F Verstraete and H. Verschelde, Phys. Rev. Lett. 90, 
097901 (2003); K. Audenaert, M.B. Plenio, and J. Eisert, Phys. 
Rev. Lett. 90, 027901 (2003); Y.C. Eldar, M. Stojnic, and 
B. Hassabi, Phys. Rev. A 69, 062318 (2004); B . Synak, K. 
Horodecki, and M. Horodecki, quant-ph/0405149 

[7] K. Audenaert, quant-ph/0402076 Proceedings Sixteenth Inter- 
national Symposium on Mathematical Theory of Networks and 
Systems (MTNS2004), Catholic University of Leuven, Bel- 
gium, 5-9 July 2004. 

[8] A.C. Doherty, PA. Parrilo, and FM. Spedalieri, Phys. Rev. 
Lett. 88, 187904 (2002); A.C. Doherty, PA. Parrilo, and FM. 
Spedalieri, Phys. Rev. A 69, 022308 (2004) 

[9] F.G.S.L. Brandao and R.O. Vianna, quant-ph/0405008 
quant-ph/0405063 quant-ph/0405096 
[10] L.M. Ioannou, B.C. Travaglione, D.C. Cheung, and A.K. Ekert, 

quant-ph/0403041 
[11] N.Z. Shor, Soviet Journal of Circuits and Systems Sciences 25, 
1 (1987). 

[12] M. Kojima and L. Tuncel, Mathematical Programming 89, 79 
(2000); A. Takeda, K. Fujisawa, Y. Fukaya, and M. Kojima, 
Journal of Global Optimization 24, 237 (2002); M. Kojima, S. 
Kim, and H. Waki, J. Ope. Res. Soc. Japan 46, 2 (2003). 

[13] J.B. Lasserre, SIAM J. Optimization 11, 796 (2001). 

[14] D. Henrion and J.B. Lasserre, IEEE Control Systems Magazine 
24, 72 (2004). 

[15] P.A. Parrilo, Structured semi-definite programs and semi- 
algebraic geometry methods in robustness and optimization 
(PhD thesis, California Institute of Technology, Pasadena, 
2000). 

[16] L. Vandenberghe and S. Boyd, Semidefinite programming. 
SIAM Review 38, 49 (1996); C. Helmberg, Semidefinite pro- 



gramming, European Journal of Operational Research 137, 461 

(2002) . 

[17] For any semi-definite program (functioning as the primal prob- 
lem), in its most general form being given by 

minimize c T x, (87) 

subject to F + £)* =1 x s F s > 0, 

one can formulate the Lagrange-dual problem, which is again a 
semi-definite problem. It is given by 

maximize -ti[F Z], (88) 

subject to tr[F 3 Z] = c s , s = 1, t, 
Z > 0. 

The key point of Lagrange duality is that any solution of the 
dual problem is a lower bound to the optimal solution of the pri- 
mal problem. This is what is referred to as weak duality. Under 
certain conditions (in particular, if there is a solution x satisfy- 
ing Fo + 2~^s=i x sF s > 0, the optimal values of the dual and 
the primal problem are identical. In this case, which is rather the 
typical one, one refers to strong duality. The idea of Lagrange 
duality is a powerful tool to formulate rigorous lower bounds to 
solutions of optimization problems. 

[18] N.S. Jones and N. Linden, quant-ph/0407 1 17 

[19] As pointed out before, for qubit systems the constraints can fur- 
ther be simplified by merely requiring 

tr[tr A3 [P«] 2 ] = (tr[P (l) ]) 2 (89) 
tr[tr Aj [P (!) ]] = tr[P W ] (90) 

forj £ /. 

[20] R.T Rockefellar, Convex analysis (Princeton University Press, 
Princeton, 1970). 

[21] This notion can be sharpened by employing the notion of weak 
membership 12211 . This can be phrased as follows: we denote for 
any convex set S C Q m and any rational 5 > with B(S, S) 
the set of all x £ Q m for which there exists ay £ S such that 

||x - y\\ 2 < 5, (91) 

and with B(S, —S) the set of all x £ S for which y € S for 
all y £ (Q m with \\x — y\\z < S. So clearly, S is a strict sub- 
set of B(S, 5), and B(S, —5) is a strict subset of 5*. The weak 
membership problem allows for two alternatives: given a ratio- 
nal element x € <Q m and a rational number S > either (i) 
assert that x £ B(S, 8), or (ii) assert that x B(S, —S). 
[22] M. Grotschel, L. Lovasz, and A. Schrijver, Geometric algo- 
rithms and combinatorical optimization (Springer, Heidelberg, 
1988). 

[23] B.M. Terhal, Phys. Lett. A 271, 319 (2000). 
[24] M. Barbieri, F. De Martini, G. Di Nepi, P. Mataloni, G.M. 
D'Ariano, and C. Macchiavello, Phys. Rev. Lett. 91, 227901 

(2003) ; M. Bourennane, M. Eibl, C. Kurtsiefer, S. Gaertner, H. 
Weinfurter, O. Guhne, P. Hyllus, D. BruB, M. Lewenstein, and 
A. Sanpera, Phys. Rev. Lett. 92, 087902 (2004). 



14 



[25] M. Curty, M. Lewenstein, and N. Liitkenhaus, Phys. Rev. Lett. 
92, 217903 (2004); M. Curty, O . Giihne, M. Lewenstein and N. 
Liitkenhaus, quant-ph/0409047 

[26] K. C hen and L A. Wu, Phys. Rev. A 69, 022312 (2004); G. 
Toth, quant-ph/0406061 

[27] M. Lewenstein, B. Rraus, J.I. Cirac, and P. Horodecki, Phys. 
Rev. A 62, 052310(2000). 

[28] It was shown in Ref. I29II that those states cannot be destilled, 
which is why they are referred to as being bound entangled. For 
more than two parties, there exist states which have a NPPT 
with respect to some splitting which are nevertheless bound en- 
tangled |U. 

[29] M. Horodecki, P. Horodecki, and R. Horodecki, Phys. Rev. Lett. 
80, 5239 (1998). 

[30] W. Diir, J.I. Cirac, and R. Tarrach, Phys. Rev. Lett. 83, 3562 
(1999). 

[31] G. Kimura, Phys. Lett. A 314, 339 (2003); M.S. Byrd and N. 

Khaneja, Phys. Rev. A 68, 062322 (2003). 
[32] P. Hyllus, C. Moura Alves, D. BruB, and Ch. Macchiavello, 

Phys. Rev. A 70, 032316 (2004). 
[33] CH. Bennett, S. Popescu, D. Rohrlich, J.A. Smolin, and A.V. 

Thapliyal, Phys. Rev. A 63, 012307 (2001). 
[34] N. Linden, S. Popescu, B. Schumacher, and M. Westmoreland, 

quant-ph/99 12039 E.F. Galvao, M.B. Plenio, and S. Virmani, 

J. Phys. A 33, 8809 (2000); S. Wu and Y. Zhang, Phys. Rev. 

A 63, 012308 (2001); A. Acfn, G. Vidal, and J.I. Cirac, Quant. 

Inf. Comp. 3, 55 (2003). 
[35] H. Barnum and N. Linden, J. Phys. A 34, 6787 (2001). 
[36] T.-C. Wei and P.M. Goldbart, Phys. Rev. A 68, 042307 (2003). 
[37] V. Coffman, J. Kundu, and W.K. Wootters, Phys. Rev. A 61, 

052306 (2000); J. Eisert and H.J. Briegel, Phys. Rev. A 64, 

022306 (2001); D.A. Meyer and N.R. Wallach, J. Math. Phys. 

43, 4273 (2002); F. Verstraete, J. Dehaene, and B. De Moor 

Phys. Rev. A 68, 012103 (2003); A.J. Scott, Phys. Rev. A 69, 



052330 (2004); M. Hein, J. Eisert, and H.J. Briegel, Phys. Rev. 
A 69, 062311 (2004). 

[38] A. Shimony, Ann. NY Acad. Sci. 755, 675 (1995). 

[39] L.-M. Duan, G. Giedke, J.I. Cirac, and P. Zoller, Phys. Rev. 
Lett. 84, 2722 (2000); J.I. Cirac, J. Eisert, G. Giedke, M. 
Lewenstein, M.B. Plenio, R.F. Werner, and M.M. Wolf, text- 
book in preparation (2004). 

[40] H.F. Hofmann and S. Takeuchi, Phys. Rev. A 68, 032103 
(2003); G. Toth, Phys. Rev. A 69, 052327 (2004). 

[41] O. Giihne, Phys. Rev. Lett. 92, 117903 (2004). 

[42] R. Horodecki, P. Horodecki, and M. Horodecki, Phys. Lett. A 
210, 377 (1996); R. Horodecki and M. Horodecki, Phys. Rev. 
A 54, 1838 (1996); K.G.H. Vollbrecht and M.M. Wolf, J. Math. 
Phys. 43, 4299 (2002); G. Adesso, F. Illuminati, and S. De 
Siena, Phys. Rev. A 68, 062318 (2003); C. Moura Alves and 
D. Jaksch, Phys. Rev. Lett. 93, 110501 (2004). 

[43] G.G Amosov, A.S. Holevo, and R.F. Werner, Problems in 
Information Transmission 36, 25 (2000); C. King and M.B. 
Ruskai, quant-ph/0401026 

[44] H.D. Sherali and W.P. Adams, SIAM Journal on Discrete Math- 
ematics 3,411 (1990). 

[45] D. Henrion and J.B. Lasserre, ACM Transactions on Math- 
ematical Software 29, 165 (2003). See also the web page 
www.laas.fr/~henrion/software/gloptipoly/gloptipoly.html. 

[46] J.F. Sturm, Optimization Methods and Software 11, 625 (1999); 
see also the documentation of the software on the web page 
fewcal.kub.nl/sturm/software/sedumi.html. 

[47] W. Diir, G. Vidal, and J.I. Cirac, Phys. Rev. A 62, 062314 
(2000). 

[48] J. Kempe, D. Bacon, D.A. Lidar, and K.B. Whaley, Phys. Rev. 

A 63, 042307 (2001). 
[49] A.C. Doherty, PA. Parrilo, and F.M. Spedalieri, 

quant-ph/0407143 



